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We consider the modification of extragalactic cosmic ray spectrum caused by cos- 



mic ray interactions with infrared background photons which are present in the ex- 
Q-r tragalactic space together with relic photons. It is assumed that cosmic ray spectrum 

at superhigh energies has extragalactic origin and is proton dominated. 
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Presence of cosmic infrared background radiation (CIB) in extragalactic space has now 
been confirmed by experiments, and its intensity is measured (see, e.g., reviews [1, 2]), 



^ ■ though the precision of these measurements is still not so high, especially if to compare with 
those for relic radiation. Nevertheless, we know enough to estimate some effects caused 



by the presence of the extragalactic infrared background. One of such effects is neutrino 
production due to high energy cosmic ray proton interactions with CIB photons, considered 
in [3-5]. Here we will focus on the other point: high energy cosmic ray spectrum modification 
■ due to interactions with CIB photons. 



2. EXTRAGALACTIC COSMIC RAY PROTON SPECTRUM 

For the calculation of cosmic ray (CR) energy spectrum at high energies, we assume it to 
be extragalactic and proton dominated from the energy E = 3 x 10 17 eV. The sources of these 
CRs are assumed to be distributed isotropically and uniformly throughout the Universe. In 
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this case it is very convenient to use cosmological transport equation written as 

MM + JL im z)n(E , z)] _ MfA =g(EtZ) . (1) 

We work in in the continuous energy loss approximation introduced in [6] and neglect possible 
proton absorbtion. In (1), n(E, z) is the number density of CR protons with a given redshift 
z, the function (3(E, z) is the change of proton energy in unit interval of z, 

P {E > Z) = lz- = —z- mZ) d-z> 

where the first term in r.h.s. is due to adiabatic energy losses (caused by cosmological ex- 
pansion), b(E, z) = —dE/dt, and is simply connected with the proton cooling rate ^(E, z). 
For the case of proton interactions with photon gas it is given by the formula 
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t-\E,z) = -b{E,z) = — 2 J de r a(e r )f(e r )e r J de- ' 
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where 7 P = E/m p is proton Lorenz factor, e r is the photon energy in the CR proton rest 
system, a(e r ) is the photoabsorbtion cross section, f(e r ) is the average relative proton energy 
loss in p7-collision (in the observer system), e t h is the threshold of the photoabsorbtion 
reaction. 

Protons lose their energy interacting with photons via jry — > nX and jry — ► pe + e~ reac- 
tions. Here, we suppose that 7 can be relic photon (they have black body spectral distribution 
with T w 2.1K) or background infrared photon. 

Thus, in our approximation the total cooling rate should be written as a sum of two com- 
ponents: cooling due to interactions with relic and infrared components of the extragalactic 
radiation background, 

t;\E, z) = t-l^E, z) + t;} nfr (E, z). 

The function g(E, z) in r.h.s. of the kinetic equation (1) describes the combined source 
of extragalactic cosmic rays. This source function can be written in the form 

g(E,z) = p(z)r ] (z)f(E)^. (2) 

Here, p(z) is the number density of local CR sources (e.g., AGNs) in the proper (physical) 
volume, p(z) = po(l + z) 3 r)(z) is the activity of each local source (the integrated number of 
produced particles per second), r)(z) = (1 + z) m r] 9(z max — z). Writing this, we assume that 



3 



the cosmological evolution of cosmic ray sources can be parametrized by power law with the 
sharp cut-off at some epoch with redshift z max (m and z max are considered as parameters of 
a model of the combined source). At last, the function f(E) in eq.(2) describes a form of the 
differential energy spectrum of the local source. For a simple power law injection spectrum 
we can write 

Eq \EqJ 

Now supposing some values for the presented parameters we can solve (1). One more thing 
needed is infrared photon number density for different redshifts, n IR (E 1 , z), to calculate 
energy losses in infrared background. As for the case of relic photons, the dependence 
of n rellc (E 1 ,z) on z is trivial: they always have black body spectrum, with temperature 
T(z)=T(0)(l + z). 

3. COSMIC INFRARED RADIATION BACKGROUND 

Extragalactic infrared background had been formed by infrared-luminous galaxies in the 
late Universe. For the calculation of CIB it is again convenient to use the cosmological 
transport equation which is analogous to that used in the previous section. The function 
which must be found is the number density of infrared photons at different cosmological 
epochs, n IR (E 1 , z). 

The resulting expression for the number density of infrared photons in extragalactic space 
is [4] 

z 

where L is the luminosity of infrared-bright galaxy emitting photons, S IR (E^,L) describes 
the spectrum of this radiation (it is called the spectral energy distribution), p(z,L) is the 
number density of infrared-bright galaxies with a given luminosity. 

From direct measurements we can learn p(0,L) and S IR (E 1 , L). Then we can assume 
some kind of the cosmological evolution to obtain p(z, L). Next step is to calculate n IR {E 1 , 0) 
making use of formula (3), and compare it with the measured CIB intensity. One of the 
results of such estimations is the well-known fact [1, 2] that we need to assume strong 
evolution of background sources in the far infrared region to explain observed CIB intensity 
(say, vl{y) ~ 30nWm~ 2 sr^ 1 at A ~ 100pm). 
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In fact, we include in our calculations only far and partially middle infrared background 
regions, not considering near infrared and optic regions. For the present calculations, we 
used the same parametrization of the function p(z, L) and same input data as described in 
[5]. But we should mention that the effect under discussion (CR spectrum modification) 
does not depend greatly on the type of infrared sources evolution, but much more on the 
CIB intensity at z = (that is, at our time). 

4. RESULTS AND CONCLUSIONS 

On Fig. la we show the proton cooling rate in infrared photon gas (together with corre- 
sponding function for the relic photon case), both curves shown for z = 0. The asymptotic 
value of tp 1 at highest energies is proportional to the total number density of photons, so 
for E > 10 20 eV cooling on infrared photons gives a small contribution, ~ 0.5%. But one can 
see from the figure that at proton energy ~ (3 4- 4) x 10 19 eV, the contribution of infrared 
photons in total cooling rate is noticeable. 

For calculation of the CR spectrum we used the following set of parameters determining 
the spectrum slope and cosmological evolution of CR sources: 7 = 2.5, m = 3.5,z max = 5. 
We assumed here that extragalactic spectrum of CRs dominates beginning from E = 3 x 
10 17 eV. The theoretical CR spectrum at z = was normalized on experimental data; from 
such a normalization we obtained the value of the product po^o entering eq.(2) for the source 
function of the CR proton kinetic equation, p ^o ~ 1-4 x 10~ 42 cm~ 3 s _1 . 

The results of two calculations of extragalactic CR proton spectrum (with and without 
taking into account energy losses on infrared photons) are shown on Fig. lb. Both curves 
have characteristic dip and cut-off features; one can see that extragalactic CR spectrum near 
GZK cut-off gets modified when taking into account interactions with infrared background. 

In this work we used concrete model parameters for CR source properties and the certain 
model for CIB, for the estimation of CR extragalactic spectrum modification at 10 19 — 10 20 eV. 
Certainly, new experimental data for CIB intensity as long as improvement of our knowledge 
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Figure 1. a) Proton cooling rates in relic and infrared photon background, 

b) Extragalactic CR proton spectrum calculated with (lower curve) and with- 
out (upper curve) taking into account energy losses in the infrared background. 

of CR sources will help to determine more precisely the magnitude of this effect. 
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